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Abstract 



We study the fragmentation of the b quark in top decay in NLO QCD, 
within the framework of perturbative fragmentation, which allows one to 
resum large logarithms ~ log(m^/m^). We show the 6-energy distribution, 
which we compare with the exact 0{as) result for a massive b quark. We 
use data from e"'"e~ machines in order to describe the 6-quark hadronization 
and make predictions for the energy spectrum of 6-flavoured hadrons in top 
decay. We also investigate the effect of NLL soft-gluon resummation in the 
initial condition of the perturbative fragmentation function on parton- and 
hadron-level energy distributions. 
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1 Introduction 



For the sake of performing accurate studies of the top-quark properties and a precise 
measurement of its mass at the present Run II of the Tevatron accelerator and, in future, 
at the LHC 0] and at the Linear Collider 0, a reliable description of the bottom-quark 
fragmentation in top decay t —>■ bW will be essential. 

As shown in ^ , the 6-fragmentation is indeed one of the sources of uncertainty in the 
measurement of the top mass at the Tevatron, as it contributes to the so-called Monte 
Carlo systematics. At the LHC, recent studies 0] have suggested that final states with 
leptons and J/ijj, with the J/ijj coming from the decay of a 6-flavoured hadron and the 
isolated lepton from the W decay, will be a promising channel to reconstruct the top 
mass. In the expected experimental error, a result of statistics and systematics, has 
been estimated to be Am^ ~ 1 GeV and the b fragmentation is the largest source of 
uncertainty, accounting for about 0.6 GeV. 

Available tools to describe the 6-quark hadronization are Monte Carlo event genera- 
tors, such as HERWIG ||^, PYTHIA or ISA JET [0, implementing respectively cluster 
P], string 1^ and independent-fragmentation [|1^ models. Monte Carlo programs de- 
scribe the initial- and final-state multiparton radiation in hadron collisions according to 



the soft and/or collinear approximation (see, for example, JTll). HERWIG and PYTHIA 
parton showers have been provided with matrix-element corrections for a few processes, 
such as top decay [Q, in order to allow hard and large-angle parton radiation. The 



analysis of |J] was in fact performed using the PYTHIA event generator. In [0, the 
HERWIG event generator was used to perform studies on the top mass reconstruction 
at the LHC, relying on the 6-quark fragmentation. The mse invariant mass distribution, 
where B is a 6-flavoured hadron and i a lepton from the W decay, was exploited in order 
to fit the top mass. 

In this paper we analyse the 6-quark fragmentation in top decay in the framework of 
perturbative fragmentation at next-to-leading order (NLO) in QCD. The factorization 
theorem ||I4| dictates that, up to power corrections ~ 0{{1/Qy), with p > 1 and Q being 
a characteristic energy scale of the process, a hadron-level cross section can be written 
as the convolution of a short- distance, perturbative cross section and long-distance, non- 
perturbative terms, corresponding to initial-state parton distribution functions and/or 
final-state fragmentation functions. For heavy-quark production, the quark mass m 
acts as a regulator for the collinear singularity and allows one to perform perturbative 
calculations. However, fixed-order event shapes or differential distributions typically 
contain terms like asloglQ"^ /rn^), where Q is, for example, the centre-of-mass energy 
or the heavy-quark transverse momentum. Such terms spoil the convergence of the 
perturbative expansion and make fixed-order calculations unreliable once Q is much 
larger than m. The method of perturbative fragmentation functions, originally proposed 
in 



15|| , allows one to resum these large logarithms. 



According to the method in |]T3[, heavy quarks are first produced at large transverse 
momentum m ^ Pr, as if they were massless, and afterwards they slow down and 
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fragment into a massive object. The perturbative fragmentation function D{^FiTt^) 
expresses the transition of a massless parton into a massive quark at the factorization 
scale hf- 

The value of D{fiF,m) at any scale fip can be obtained by solving the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations |jl6|,|T^, once its initial value at 



a scale /iof is assigned. The universality of the initial condition and, in general, of 
the perturbative fragmentation function, already suggested in in the framework 
of e~^e~ annihilation, has been recently proved in a completely process-independent 
way [jl8[. As discussed e.g. in |19| for heavy-quark production at hadron colliders, the 



perturbative fragmentation formalism yields a weaker dependence of phenomenological 
observables on the renormalization/factorization scales and on the chosen set of parton 
distribution functions. Furthermore, the analysis of [^] fully resums the leading (LL) 
and next-to-leading logarithms (NLL) which are associated with the emission of soft 
gluons and appear in the initial condition of the perturbative fragmentation function 
(process independent) and in the parton-level differential massless cross section (process 
dependent) of e~^e~ annihilation. 

Finally, in order to describe the non-perturbative fragmentation of a parton into a 
hadron, several phenomenological models have been proposed pO|j2T[| , besides the ones 
which are implemented in Monte Carlo event generators. Non-perturbative fragmen- 
tation functions contain parameters which need to be fitted to the experimental data. 
Since the hadronization mechanism is universal and independent of the perturbative pro- 
cess which produces the heavy quark, one can exploit the existing data on e+e~ —>■ bb 
events to fit such models and describe the 6-quark non-perturbative fragmentation in 
other processes, such as top decay. 

Perturbative fragmentation functions and non-perturbative hadronization models 
have been extensively applied to study the physics of c- and/or 6-fiavoured hadrons 



produced in e^e annihilation p2| - p5| , hadron collisions [p^p6| . Deep Inelastic Scattering 
p4| and 7p collisions 



In this work, we apply this method to the 6-fragmentation in top decay. In Section 
2 we review the method of perturbative fragmentation functions and apply it to predict 
the 6-quark energy spectrum in top decay. In Section 3 we analyse the non-perturbative 
fragmentation of the b quark and show energy distributions of 6-hadrons in top decay, 
making use of fits to LEP and SLD data to parametrize the hadronization models. 
In Section 4 we summarize the main results of our analysis and make comments on 
possible developments of our study. In Appendices A and B we show details of our 
calculation and, comparing results for massless and massive b quarks, we check that 
our computation is consistent with the initial condition of heavy-quark perturbative 
fragmentation functions. 
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2 Perturbative fragmentation and parton-level re- 
sults 



We wish to study 6-quark production in top decay within the framework of perturbative 
fragmentation functions. We consider the decay of an on-shell top quark at next-to- 
leading order in as Q 

t{q) ^ h{p,)W{pw) {9{Pg)) (1) 
and define the 6-quark scaled energy fraction 

Xe = 2~- (2) 

mi 

Neglecting the b mass, we have < xe < ^ — w, w being w = rriy^/mf. Throughout 
this paper, we shall make use of the normalized b energy fraction: 

Xb = — , <Xb <l. (3) 
1 — w 

Following , the differential width for the production of a massive b quark in top decay 



can be expressed via the following convolution: 

1 dr 1 dz dTi f^b 



Tq dxi 



ixb,mt,mw,mb) = f —-r^{z,mt,mw,f^,IJ'F)Di(—,iJF,mb\ (4) 
To i Jxt z dz \ z J 



where Fq is the width of the Born process t — > bW ^ dVi/dz is the differential width for the 
production of a massless parton i in top decay with energy fraction z, Di{xb/ z, fiF,mb) 
is the perturbative fragmentation function for a parton i to fragment into a massive b 
quark. In Eq. (^) fi and fiF are the renormalization and factorization scales respectively, 
the former associated with the renormalization of the strong coupling constant. In 
principle, one can use two different values for the factorization and renormalization 
scales; however, a choice often made consists of setting fi = fiE and we shall adopt this 
convention for most of the results which we shall show. 

The approach of and the factorization on the right-hand side of Eq. (^) are 



rigorously valid if one can neglect terms behaving like {m/Qy, where p > 1, m is the 
mass of the fragmenting heavy quark and Q is the hard scale of the process. This is 
indeed our case since the scale of top decay is set by its mass and mb/rrit — (9(10~^). 

The definitions of dVi/dz and Di{xb/ z, fiFymb) are not unique, but they depend on 
the scheme which is used to subtract the coUinear singularities which appear in the 



massless differential width dVi/dz. In [|15[, the MS factorization scheme is chosen and 
we shall stick to it hereinafter. 



^Our assumption Bit hW) = 1 is consistent with recent measurements of the CDF Collaboration 
of the ratio R = B{t Wb)/B{t Wq), where g is a d, s or 6 quark, and the subsequent extraction 
of the Cabibbo-Kobayashi-Maskawa matrix element Vtb [p3 ■ 
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Since we have been assuming B(t ^ bW) = 1 and the probabihty to produce a b 
quark via the sphtting of a secondary gluon is neghgible, we shall safely limit ourselves 
to considering the perturbative fragmentation of a massless b into a massive b and, on 
the right-hand side of Eq. (^), we shall have only the i = b contribution. We shall 
then need to evaluate the differential width dTb/dz for the production of a massless b 
quark in top decay, in dimensional regularization, and subtract the collinear singularity 
according to the MS prescription. We obtain the following MS coefficient function: 



l_dl\ 
Tq dz 



MS 



(5(1 - ^) + — Ai{z), 



2tx 



(5) 



with 

AAz) 



+ 



Cf 
l + z' 

2(1 + z^ 



1 + z' 



-6{l-z) 



l-z)+ ■ 2 
4 log [(1 - w)z\ 
1 



m^ ^1 



log^ + 2 



w 



1 



l + 2w 



log(l - z) 



+ 6{l-z 
2w 



l + 2w 

Az 

■ logz 

27r2 



-21og(l-w) 

w{l-w){l- zf 
'l + 2w){l- {l-w)z) 



4Li2(l - w) + 21og(l -w)\ogw — + 



l + 2w 



log(l — w) 



1 — w 



log w + 



3w 



l + 2w 



-9 



(6) 



where 



Li2(x) 



^ df 
t 



(7) 



is the Spence function. In Appendix A we shall give more details on the derivation of 
Eq. (H) and present results for the differential width dV/dxh once the 6-quark mass is 
fully taken into account. 

In order to be consistent at NLO, in Eq. @ «5(/i) is to be the strong coupling 
constant at NLO as well: 



with 6o and bi given by 



6olog(Ai7A2 



&ilog[log(/iVA^)] 
62log(/iVA2) 



33 - 2n 



f 



127r 



153 - 19n 



f 



247r2 



(9) 



A being the typical QCD scale and = 5 the number of flavours. 

We note that the coefficient function contains a term where the strong coupling 
constant multiplies the logarithm log{m^ / fj?p) . For our calculation to be reliable, we shall 
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have to require such a logarithm not to be too large, which implies that the factorization 
scale hf will have to be chosen of the order of rrit. 

In [15 1, considering heavy-quark production in e+e" annihilation and comparing 
the massive and massless differential cross sections in the MS factorization scheme, 
the authors have obtained the NLO initial conditions at a scale /iof for heavy-quark 
perturbative fragmentation functions. For the 6 — > 6 transition, it has been found [|: 



27r 



1- Xb 



flog^-21og(l 



Xb) - 1 



(10) 



with Ci? = 4/3. In order to avoid large logarithms in Eq. (|T0D, the scale /xof is to be 
taken of the order of the b mass. The universality of the initial condition ([To|) has been 
lately proved in 



1^ 



and one can therefore exploit it to predict the b fragmentation in 
top decay as well. 

For the sake of completeness, in Appendix A we shall show the result for the differ- 
ential width dr/dxb once we keep the b mass only in contributions ~ log(m^/m^) and 
neglect terms proportional to powers of the ratio mb/rrit. Comparing the result with 
the massless rate (^, we shall be able to reproduce the initial condition of the 6-quark 
fragmentation function, which will be a consistency check of our calculation and, at the 
same time, a confirmation of the validity of Eq. ([T0|) in our context as well. 

Assigned the initial condition (0), the value of the perturbative fragmentation func- 
tion at any other scale fip can be obtained by solving the DGLAP evolution equations 

nini: 



d 



d log jjp 



Di{xb,HF,rnb) = J2f ~Pio Oisdj-F)) Dj{z, HF^rnb), 

• Jxt, Z \ Z J 



:iii 



where 



Pij\^b) are the Altarelli-Parisi splitting functions ||T6[, and the higher-order terms 



Pij{xb,as{fiF)) 



27T 



(Xb) + 



27r 



P^\xb) + 0{a 



(12) 



Pl:l\xb) can be found in |9|,|0| 



As shown in |]T5[, solving the DGLAP equations for the evolution from a scale /^qf 
to /ii? allows one to resum potentially-large logarithms ~ cx^^f) ^'^^{i^^f I iAf) ■ Assuming 
/ioF — f^h and ^F — ^t, and considering the sphtting functions (|T^) at 0{as), one 
resums the leading logarithms 



Q;g(mt, 



Eq. (0) leads to the inclusion of next-to-leading logarithms 



\og"'{mf/ml). Accounting for C(a|) terms in 



a 



n+l. 



\og^ {mf /ml) as 



well. In this paper, we shall always assume that the 6-quark perturbative fragmentation 
function evolves with NLL accuracy. 



^For the fragmentation of a gluon, light quark or b quark into a b quark, see [h5| 
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The DGLAP equations get highly simphfied in Melhn space; the solution reads 



Di^N{lJ'F,rnb) 



+ 



A,Af(/ioF, rrib) exp l^P^h 
1 



p(i) _ 27r6i (0) 



(13) 



with Di^^ifJ'F, f^b) being the Melhn transform of the x-space perturbative fragmentation 
function 

Di^N{lj,F,mb) = dx Xb'^Di{xb, lJ,F,rnb) (14) 

JO 

and the variable t defined as 



27r6f 



log 



OF) 



(15) 



Throughout our analysis, i = b in Eqs. (0-|T^) and Db,N{l^oF,fnb) is the iV-space trans- 
form of Eq. (|T0|). Expressions for Db^N{fJ'OF,fnb) and for the Mellin transforms of the 
NLO splitting functions in Eq. (0) P^y-* and P^^^ can be found in [|I5|. We shall have 
to compute the iV-space transform of the MS coefficient function (|^) f\ and multiply it 
by Eq. (O) in order to get the Mellin transform of Eq. (0): 



TN{mt, mw, mb) = f Ar(mi, mw, yU, fiF)Db,N{fJ'F, rrib) 



(16) 



with 



r 



N[mt,mw,mb) = — 



dr 



{Xb,mt,mw,rnb) 



(17) 



1 

To ^0 " dxb 

The 6-quark energy distribution in x-space will finally be obtained by inverting the 
A^-space result (|T^) numerically. 

Furthermore, in Eq. (|IUp the coefficient multiplying the strong coupling constant 
contains terms behaving ~ 1/(1 — Xb)+ or ~ [log(l — Xb)/{1 — Xb)]+ once Xb — > 1, 
which corresponds to behaviours ~ log or ~ log^ in moment space, for large A^. 
The limit Xb ^ I {N ^ oo) corresponds to soft-gluon radiation in top decay. Soft LL 
~ ag(/io) log"""*"^ A^ and NLL ~ «s(Aio) log" A^ contributions in the initial condition of the 
perturbative fragmentation function have been resummed in [|18[. Due to the process 



independence of the heavy-quark perturbative fragmentation function, we can exploit 
the result of []T8[ , which we do not report here for the sake of brevity, to resum LL and 



NLL terms in Eq. (13) in top decay as well. 

We wish to present results for the normalized 6-quark energy distribution in top decay 
using the technique just described. We normalize our plots to the total NLO width F, 
obtained neglecting powers ~ (mb/mj)^, whose expression can be found in In 



^In Appendix B we shall present results for the A^-space counterpart of Eq. (^) . 
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fact, the factorization on the right-hand side of Eq. the DGLAP evolution and the 
resummation of soft logarithms in the initial condition of the perturbative fragmentation 
function with NLL accuracy do not affect the total NLO normalization^. 

It can be observed that, as long as one neglects interference between top production 
and decay, one has: 

1 da 1 dr , , 

a dxb 1 dxb 

where {l/a)da/dxb is the normalized differential cross section for the production of 
a b quark with energy fraction Xb via top quarks, independently of the production 
mechanism. Our results will then be applicable to pp (Tevatron), pp (LHC) or e~^e~ 
(Linear Collider) collisions. For our numerical study, we shall assume = 175 GeV, 
mw = 80 GeV, mb = 5 GeV and A = 200 MeV. 

We show our results for the Xb spectrum in Fig. |I[ We plot the Xb distribution 
according to the perturbative fragmentation approach, with and without NLL soft- 
gluon resummation in the initial condition of the perturbative fragmentation function. 
For the sake of comparison, we also show the exact 0{as) result for a massive b quark, 
whose analytical expression will be given in Appendix A. We set fi = fip = and 
1^0 = fJ-oF = rUb. 

We note that the use of perturbative fragmentation functions has a stronger im- 
pact on the Xb distribution. The fixed-order result lies well below the perturbative 
fragmentation results for about 0.1 <Xb< 0.9 and diverges once Xb 1, due to a be- 
haviour ~ 1/(1 —Xb)+- Moreover, the full inclusion of powers of mb/rrit has a negligible 
effect on the Xb spectrum; the dot-dashed and dotted lines in Fig. |l| are in fact al- 
most indistinguishable. As for the perturbative fragmentation results, the distribution 
with no soft-gluon resummation shows a very sharp peak, though finite, once Xb ap- 
proaches unity. This behaviour is smoothed out after we resum the soft NLL logarithms 
appearing in the initial condition of the perturbative fragmentation function, as the 
6-energy spectrum gets softer and shows the so-called 'Sudakov peak'. Both perturba- 
tive fragmentation distributions become negative for a;;, ^ and Xb —>■ 1, which is a 
known result, already found for heavy-quark production in e+e~ annihilation [p!^ ,p5[]. 
For Xb — >■ 0, the coefficient function (|^) contains large logarithms ~ ^^(/i) logx;, which 
have not been resummed yet. Likewise, in the soft limit Xb ^ 1, Eq. (^ contains contri- 
butions ~ Ois{n)/ {1 — Xb)+ and ~ ^^(/i) [log(l — X6)/(l — Since asijnb) — 2as{mt), 
for fi = nit and /iq = mb such terms are smaller than the similar ones which appear 
in the initial condition of the perturbative fragmentation function ([T0|) , but nonetheless 
they would need to be resummed. As stated in [113, once Xb gets closer to unity, non- 



perturbative contributions also become important and should be taken into account. 
The region of reliability of the perturbative calculation at large Xb may be related to 
the Landau pole in the expression for the strong coupling constant (§), and estimated 



^As far as soft-gluon resummation is concerned, matching the resummed initial condition to the exact 
0{as) result (foh (see Eq. (76) in Ref . [jisl ) , guarantees that the total NLO width is left unchanged. 
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Figure 1: 6-quark energy distribution in top decay according to the perturbative frag- 
mentation approach, with (sohd hne) and without (dashes) NLL soft-gluon resummation 
in the initial condition of the perturbative fragmentation function, and according to the 
exact NLO calculation, with (dot-dashes) and without (dots) inclusion of powers of 
mb/nit. In the inset figure, we show the same curves on a logarithmic scale. 



to be Xf, < 1 — A/mf,. 

Figs. 1^ and ^ show the dependence of the perturbative fragmentation Xf, distribu- 
tions on the factorization scales /iof and fip, with (Fig. ^ or without (Fig. |^) soft-gluon 
resummation. We observe in Fig. ^ (a) that the dependence on the initial scale /iof is 
small when we resum soft logarithms in the initial condition of the 6-quark perturbative 
fragmentation function. Little impact is only visible in the neighbourhood of the Su- 
dakov peak and, in particular, the distributions obtained for fiQp = nih and fi^p = 2mf, 
are indistinguishable from each other. On the contrary. Fig. |^ (a) shows that the Xf, 
spectrum has a remarkable dependence on //qf once soft logarithms are not resummed. 
Comparing Fig. ^ (b) and Fig. ^ (b), we note a quite similar dependence on the choice 
oi lip- In fact, /xof rather than yUi? is the scale which enters the expression of the initial 
condition Db{xb, fioF, mh). For ^ip approaching /xof, the distribution with no soft resum- 
mation gets closer to the fixed-order, unevolved one shown in Fig. |l]. We expect that 
once one includes soft resummation in the MS coefficient function (|), which contains 
the factorization scale /ii?, the dependence of the 6-energy spectrum on ^p will become 
weaker as well, as found in |jl8| for the purpose of e^e~ annihilation. As a whole, we can 
state that resumming soft logarithms yields a reduction of the theoretical uncertainty, as 
the dependence on factorization scales is indeed an estimation of effects of higher-order 
contributions which we have been neglecting. 
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Figure 2: (a): Xb spectrum ior = nit and /iof = f^bf^ (dots), /ioF = f^b (solid) 
and /ioF = 2mfc (dashes); (b): fioF = f^b and fip = rnt/2 (dots), fip = f^t (solid) and 
= (dashes). The renormalization scales are kept at fi = rrit and fiQ = rrih. All 
curves include NLL soft-gluon resummation in the initial condition of the perturbative 
fragmentation function. 
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3 Non-perturbative fragmentation and hadron-level 
results 



In this Section we shall present results for the energy distribution of 6-flavoured hadrons 
in top decay. We consider the transition b ^ B, where B is either a meson or a baryon 
containing a b quark and define the normalized 6-hadron energy fraction xb similarly to 
the parton-level one in Eq. (^. 

The energy distribution of a hadron B can be expressed as the convolution of the 
parton-level spectrum with the non-perturbative fragmentation function (xB)'■ 



^^ {xB,rnt,mw,mb) = ^ [ — ^{z,mt,mw,rnb)D!^(—\ (19) 

i Jxb z dz \ z J 



To dxB 

In Eq. (p!9|), {l/TQ)dT / dz is the parton-level differential width for Xi, = z and 
D^p{xb/z) is the non-perturbative fragmentation function describing the hadronization 
b ^ B, which is process independent. We can therefore exploit data from e~^e~ bb 
processes to predict the 6-quark hadronization in top decay. 

Several models have been proposed to describe the non-perturbative transition from 
a quark- to a hadron-state. One of the most-commonly used p2| , |2^ , |26| . |27[| consists of a 
simple power functional form: 

B{x,y) being the Euler Beta function. 

The model of Kartvelishvili et al. is still a power law, but with just one free 
parameter 6: 

^-(-^*) = (rTiy(2T^)'^ '''' 

We expect that if we use the non-perturbative fragmentation function in (^), which 
has two fittable parameters, we shall be able to get better agreement to the data than 
when using (0). However, in our analysis we shall try to tune the model (^) as well, 
in order to investigate how good it is at reproducing the data and how it compares to 
the other models considered. 

Finally, the Peterson model |PT[ describes the transition of a heavy quark into a 



heavy hadron according to the following non-perturbative fragmentation function: 

For an explicit expression of the normalization factor A and of the Mellin transform of 
Eq. (P^, see The A^-space transforms of Eqs. pOj ) and (pi]) are straightforward. 
The parameters a and (3 in Eq. (pOD, 6 in Eq. (^) and e in Eq. (p2D will have to be 
obtained from fits to experimental data. 
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To predict the 6-quark hadronization properties, we shall use LEP data from the 
ALEPH Collaboration |Q and SLD data [0 for e"'"e~ collisions at the Z pole, i.e. 



a/s = 91.2 GeV. Both data sets refer to weakly-decaying 6-hadrons ^ however, while 
the ALEPH x^-data just accounts for B mesons, the SLD set also considers 6-flavoured 



baryons, mainly the A^. In I^Sf, it is stated that the mean values Ixb) of ALEPH 
and SLD are consistent with each other, within the range of systematic and statistical 
errors. However, no complete statistical analysis aiming at checking the consistency 
of the xb distributions has been done yet; therefore some difference is to be expected 
when comparing fits of non-perturbative models to data actually referring to different 
6-hadron samples, as we do have. 

In order for the results of our fits to be applicable to the 6-hadronization in top decay, 
we shall have to describe the perturbative process e^e~ bb{g) in the same framework 
as we did for t —>■ bW{g) when we do the fits. In fact, the factorization on the right- 
hand side of Eq. (|l^) and the splitting between perturbative and non-perturbative part 
is somewhat arbitrary and the parametrization of the non-perturbative model indeed 
depends on the approach which is used to describe the perturbative, parton-level process 
and on the values which are chosen for quantities like A, irib and for the renormalization 
and factorization scales. We shall therefore use MS coefficient functions [^] for e~^e~ 
annihilation into massless quarks and convolute them with the perturbative 6-quark 
fragmentation function evolved to NLL accuracy according to the DGLAP equations, 
possibly including soft-gluon NLL resummation in the initial condition of the pertur- 
bative fragmentation function. We shall then obtain a parton-level differential cross 
section 1 / ao^dcr / dxb) , equivalent to Eq. (H), which we shall have to convolute with the 
non-perturbative fragmentation function of the hadronization model considered. 

The experimental analyses |33,34| use the PYTHIA Monte Carlo event generator 
to simulate e+e~ — > bb processes and the subsequent parton showering; as a result, their 
framework is pretty different from the one which we have been using for top decay. We 
cannot therefore simply use the parametrizations of the hadronization models as they 
are reported in p3| , |3^ . 

An extensive phenomenological study of fragmentation in e~^e~ processes, with more 
details on fits to LEP and SLD data is currently under way [^. For the purposes of 
our paper, we just tune the hadronization models to the data sets for a particular choice 
of the quantities which enter the perturbative calculation. 

We point out that, although we convolute our perturbative result with a non- 
perturbative, smooth function, problems are still to be expected once — > 0, 1. The 
region — > will not be reliably described since the perturbative calculation itself 
includes unresummed ~ logXf, terms. For — > 1 one should in principle correctly 
account for all the missing, non-perturbative power corrections. Resumming a class of 



^By weakly-decaying 6-hadrons we mean hadrons containing a b quark which decay according to a 
weak transition. For example, in the decay chain B* Bj {^D^*Hv) 7, the B meson rather than 
the B* is the weakly-decaying 6-hadron whose energy fraction is experimentally measurable. 
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A T PPTT 


OVjU 


(y. 


u.o 1 HZ u. 1 


1 SS -1- n zL9 




13.21 ± 1.62 


27.04 ±4.02 


x2(«,/3)/dof 


2.62/14 


11.12/16 


5 


20.39 ±0.77 


18.80 ±0.60 


X^(5)/dof 


17.27/15 


17.46/17 


e 


(1.12 ±0.16) X 10-3 


(1.17±0.10) X 10-3 




22.96/15 


130.80/17 



Table 1: Results of fits to e^e" hh data, using NLO coefficient functions, NLL DGLAP 
evolution and NLL soft-gluon resummation in the initial condition of the perturbative 
fragmentation function. We set A = 200 MeV, /ioF = /^o = "^fe = 5 GeV and ftp = = 
y/s = 91.2 GeV. a and P are the parameters in the power law (pOD, 6 refers to (|2TD, e 

to dn). 

perturbative soft logarithms and using a specific functional form with few fittable pa- 
rameters to describe the non-perturbative fragmentation is not sufficient to be able to 
include all such terms. 

In order to perform trustworthy fits to the e^e~ data and acceptable predictions 
for the 6-hadron spectrum in top decay, we shall limit ourselves to analysing data not 
too close to the critical points = 0, 1. In particular, we shall consider ALEPH data 
within the range 0.18 <xb< 0.94 and SLD data for 0.18 < a;^ < 0.90, which implies that 
our predictions for top decay will have to be considered in the same xb ranges. When 
doing the fits, we account for both statistical and systematic errors on the data. 

In Tables ^ and ^ we report the parameters which correspond to our best fits to 
the data, along with the P^r degree of freedom. We also investigate the impact of 
NLL soft-gluon resummation in the initial condition of the perturbative fragmentation 
function. Standard deviations for best-fit parameters are included as well. 

From Table we learn that the use of the power law (|20D which has two tunable 
parameters leads to excellent fits to both ALEPH and SLD data, but the values of a and 
(3 which minimize the ^'^'^ affected by fairly large errors. The model of Kartvelishvili 
et al. is good at fitting in with the ALEPH and SLD data as well. The Peterson model 
is marginally consistent with the ALEPH data and unable to reproduce the SLD data. 
Although we are comparing data samples with different 6-hadron contents, we observe 
that the best-fit values of e and 6 obtained for ALEPH and SLD are compatible within 
one and two standard deviations respectively. A bigger difference between ALEPH and 
SLD is found once we try to fit the power law (|20|) . 

Comparing Tables and we observe that the implementation of soft-gluon re- 
summation in the perturbative fragmentation function results in statistically-different 
values of the best-fit parameters. With no soft resummation, using power laws with one 
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A T PPTT 


OLjU 


(y. 


u.uu HZ u. 10 


9 nc; -1- n 98 
z^.Vjo hi u.^o 




12.39 ± 1.04 


22.10 ±2.13 


x2(«,/3)/dof 


7.12/14 


40.23/16 


5 


14.97 ±0.44 


14.57 ±0.37 


X^(5)/dof 


13.30/15 


58.63/17 


e 


(2.87 ±0.21) X 10-3 


(2.33 ±0.19) X 10-3 




52.76/15 


275.69/17 



Table 2: As in Table 0, but without resumming soft logarithms. 

or two fittable parameters still yields very-good fits to the ALEPH data, while none of 
the models considered is able to describe consistently the SLD data and the Peterson 
non-perturbative fragmentation fails in reproducing the ALEPH data as well. 

We wish now to predict the 6-hadron spectra in top decay, using models and parame- 
ters which give reliable fits to the e^e~ data. In order to account for the uncertainties on 
the parameters in the non-perturbative fragmentation functions, as shown in Tables |I] 
and 0, for each hadronization model we shall plot two curves which delimit a set of 
predictions at one-standard-deviation confidence level |^. 

In Figs. ^ and ^ we show the xb distribution using all three considered models 
fitted to the ALEPH data (Fig. D and the power forms (|0D and (|21]) fitted to SLD 
(Fig. I), with all curves including NLL soft resummation in the initial condition of 
the perturbative fragmentation function. We note that different hadronization models 
yield statistically-different predictions for 6-hadron spectra in top decay, within the 
accuracy of one standard deviation. However, one can show that the xb distributions 
according to the models (|20|) and (0) fitted to the SLD data are consistent within two 
standard deviations. This result can be related to the use of similar functional forms 
and to the large errors on a and 13. The prediction obtained fitting the Peterson non- 
perturbative fragmentation function to the ALEPH data looks pretty different from the 
others, especially at small and middle values of a;^, where the predicted errors are pretty 
small for all the considered models. Moreover, the Peterson distributions are peaked at 
slightly-larger values of xb- 

In Fig. I we compare the ALEPH and SLD predictions according to the power law 
of Eq. ( [20| ) since, as shown in Table this is the only hadronization model where the 
fitted parameters are statistically different. In fact, the overall shapes of the ALEPH- 
and SLD-based distributions lead to different predictions within one-standard-deviation 
accuracy, expecially for > 0.7, with the SLD ones being peaked at smaller xb values. 

^We point out that the correlations between the errors on a and (5 of the non-perturbative fragmen- 
tation function ( po| ) are correctly taken into account throughout our analysis and in the plots which we 
show. 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 4: xb spectrum in top decay, with the hadronization modelled according to a 
power law (solid lines), the Kartvelishvili et al. (dashes) and the Peterson (dots) model, 
with the relevant parameters fitted to the ALEPH data. The plotted curves are the 
edges of bands at one-standard-deviation confidence level. NLL soft-gluon resummation 
in the initial condition of the perturbative fragmentation function is included. We set 
jjLp = = rrit and yUoF = l^o = ^b- 

This result can be checked to be true even at higher confidence level and, as antici- 
pated, can be associated with the different 6-hadron samples which the two experiments 
reconstructed. For the sake of comparison, we also show the ALEPH-based prediction 
without resumming soft logarithms, as Table ^ reports a small even in this case. 
We find that, although the perturbative predictions look rather different according to 
whether such contributions are resummed or not (see Fig. |l|), the hadron- level results 
agree for In fact, the convolution with the non-perturbative fragmentation 

function smears the sharp peak shown by the unresummed result and the parameters a 
and P are accordingly set by the fit to the e'^e~ data in such a way that the two xb- 
predictions do not differ too much from each other. Statistically-significant differences 
are nonetheless found for xb < 0.5, where the predicted bands get narrower. 

It is finally interesting to gauge the overall impact of the non-perturbative fragmen- 
tation by comparing our parton-level result with one of our hadron-level predictions, 
in particular the one obtained using the power law with two parameters fitted to the 
ALEPH data, as it corresponds to our smallest x^- From Fig. 0, we learn that the 
hadronization effects are remarkable and the xb distribution is considerably softened 
with respect to the Xb one. This means that, even when one uses a refined perturba- 
tive approach, with DGLAP evolution and soft resummation in the initial condition of 
the perturbative fragmentation function up to NLL accuracy, the role played by the 
non-perturbative input and hence by the e~^e~ experimental data will still be crucial to 
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0.0 0.2 0.4 0.6 0.8 1.0 

Xb 

Figure 5: As in Fig. ^, but fitting the hadronization- model parameters to the SLD data. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 6: As in Fig. ^, using the model of Eq. (pOD and fitting its parameters to ALEPH 
(solid lines) and SLD (dashes), including NLL soft resummation in the initial condition 
of the perturbative fragmentation function. We also show results from the fit to ALEPH, 
but with no soft resummation (dots). 
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Figure 7: Comparison of hadron- {j = B, solid) and parton-level {j = b, dashes) results 
using the hadronization model of Eq. (pOD, with a and (3 being the central values in the 
fit to the ALEPH data, and including NLL soft resummation in the initial condition of 
the perturbative fragmentation function. 

perform any prediction on 6-fragmentation in top decay. 

4 Conclusions 

We discussed the 6-quark fragmentation in top decay in NLO QCD using the method of 
perturbative fragmentation, which resums large logarithms ~ log(m^/m^) which multi- 
ply the strong coupling constant in the fixed-order massive calculation. 

We computed the NLO differential width of top decay with respect to the 6-quark en- 
ergy fraction Xb for a massless b, fully including 6-mass effects and for a massive b quark, 
but neglecting contributions proportional to powers of the ratio nib/ nit. We determined 
the MS-subtracted coefficient function and checked that our result is consistent with 
the known expression of the initial condition for heavy-quark perturbative fragmenta- 
tion functions. We convoluted our MS coefficient function with the process-independent, 
perturbative fragmentation function for a massive b quark, evolved up to NLL accuracy 
using the DGLAP equations. We showed results for the b energy-fraction distribution 
in top decay, which we compared to the fixed-order results for a massive b quark. We 
found that the use of the perturbative fragmentation approach has a remarkable effect 
on the parton-level distribution which is smoothed out with respect to the 0{as) one, 
which gets arbitrarily large once Xh approaches unity. We also investigated the impact 
of the resummation of process-independent next-to-leading logarithms, which appear in 
the initial condition of the perturbative fragmentation function Db{xb, fioF,^b) and are 
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associated with soft-gluon radiation. We found that it softens the Xb distribution and 
makes the dependence on the scale /iof weaker. 

We then studied the energy distribution of 6-hadrons in top decay and fitted some 
hadronization models to e~^e~ data. We used ALEPH data on 6-fiavoured mesons and 
SLD data on 6-fiavoured baryons and mesons. In order to perform such fits we described 
the perturbative process e"'"e~ — > bb{g) as we did for 6-quark production in top decay. 
Throughout our analysis, we discarded data points where the hadron-level energy frac- 
tion is close to Xb — 0,1, since our approach is unreliable in the neighbourhood of these 
points. 

We found that, within our perturbative framework, models which describe the non- 
perturbative fragmentation according to power laws with one and, in particular, with 
two fittable parameters lead to good descriptions of ALEPH and SLD data. The Peter- 
son model is marginally consistent with the ALEPH results and unable to describe the 
SLD data. We also found that, although the ALEPH and SLD data refer to different 
6-hadron samples, the best-fit parameters obtained using the Kartvelishvili or the Pe- 
terson model are statistically-consistent within the error ranges. The implementation 
of NLL soft-gluon resummation in the initial condition of the perturbative fragmenta- 
tion function yields pretty different fits to the e'^e~ data. An unresummed perturbative 
calculation still leads to good fits to the ALEPH data once we use power laws to model 
the hadronization, while it is unable to reliably describe the SLD data. 

We then showed distributions of the energy fraction a;^ of 6-fiavoured hadrons in top 
decay, using only models which give rehable descriptions of the e+e" data. We found 
that the models fitted to the ALEPH data yield statistically-distinguishable predictions 
for the 6-meson spectrum in top decay. In particular, the Peterson-model distributions 
lie quite far from the others and are peaked at slightly-larger xb values. If we model 
the hadronization using power laws with one or two parameters, fitted to the SLD data 
sample, the predictions for the xb spectrum in top decay are compatible within two 
standard deviations, a result which is due to the large uncertainties on the best-fit 
values of a and j3. Wc also compared results obtained using the power law with two 
parameters, but fitted to the ALEPH or SLD data, and obtained predictions which are 
statistically different. 

We investigated the impact of NLL soft-gluon resummation in the initial condition 
of the perturbative fragmentation function on hadron-level distributions using ALEPH- 
based fits. We found a significant impact on the xb spectrum only at relatively-small 
Xb, while predictions with or without soft resummation are indistinguishable at large 
Xb values. This result stresses the importance of the non-perturbative input and conse- 
quently of the use of the e+e~ results to perform reliable predictions in top decay. This 
can be learned from direct comparison between parton- and hadron-level distributions 
shown throughout the paper. 

It will be now very interesting to use the present approach to perform predictions 
of other observables relying on the 6-fragmentation in top decay, such as the invariant- 
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mass distributions used in P,p!3| to fit the top mass value. It is clearly mandatory to 
compare the results obtained in the framework of perturbative fragmentation functions 
with the ones of Monte Carlo event generators, taking particular care about the induced 
uncertainty on the top mass measurement. However, for such a comparison to be reliable, 
Monte Carlo programs will have to be tuned to fit with the experimental data. This is 
in progress. 



Finally, we plan to extend the method developed in [|l^ in the framework of e~^e~ 
processes to resum with NLL accuracy also the process-dependent soft logarithms ~ 
a'g{fi)\og"'~^^ N and ~ ag{^)\og^ N appearing in the Mellin transform of the MS co- 
efficient function (|^). This further step will allow one to include in the perturbative 
calculation all soft logarithms, besides the process-independent ones which are present 
in the initial condition of the perturbative fragmentation function and have been cor- 
rectly accounted for throughout this paper. It will be interesting to analyse the impact 
of soft-gluon resummation in the MS coefficient function on the energy spectrum of b 
quarks and 6-flavoured hadrons. This is in progress as well. 
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A Details of the calculation 

We wish to give some details on the calculation of the top-decay differential width at 
0{as) for massless and massive b quarks. In particular, the comparison of the two re- 
sults will allow us to obtain the initial condition (|10D of the perturbative fragmentation 
function. We adopt the on-shell mass-renormalization scheme and use dimensional reg- 
ularization to regulate the ultraviolet and soft singularities and, in the massless case, the 
coUinear divergence as well. We define the parameter e, which is related to the number 
d of dimensions via d = A — 2e. 

In the massless case, the differential width will cointain a pole ~ 1/e, due to the 
collinear singularity, which disappears only in the total NLO width. This requires that, 
in order to get the correct finite term in the normalized {l/ro)drb/dxb, with Xb defined 
in Eq. (|^), the Born width Fq will have to be evaluated in dimensional regularization at 
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0(e). We find: 



J- — . ^ ■ 9 ^ i J- + e 



16 sin 6w 



w 



-7ij + log47r - 2 log(l -w) + 2 



l + w 



l + 2w 



(23) 

where a is the electromagnetic couphng constant, 7^; = 0.577216 ... is the Euler con- 
stant, Ow is the Weinberg angle and w is the ratio w = m'yy/'mf, already introduced in 
Section 2. We thus obtain: 



1 dTb 
To dxb 



i + 7ij -log47r) +Ai{xb) 



(24) 



with Ai{xb) defined in Eq. (|^). We note that Ai{xb) depends on the scale hf, remnant 
of the regularization procedure, which disappears only in the total width. We also 
checked that the integral of Eq. (^) agrees with the result of [^, where the 0{as) 
corrections to the top width have been evaluated in the approximation of a massless b 
quark. In order to get the MS coefficient function (|^), by definiteness, we shall have 
to subtract from Eq. (p^ the 0{as) term multiplying the characteristic MS constant 
(-l/e + 7£;-log47r). 

We wish to derive the NLO differential width with the full inclusion of the b mass 
mft. For this purpose, it is convenient to define the following quantities: 



P 

Q 

Go 



mi 

h^\^b-w) 
Vb 



l + b-2w + 



[I -by 



w 



xl — P'^ — arcth 



Xb 



where Xb is the normalized 6-quark energy fraction 

Xe 



Xb 



2s 



/3 < X, < 1, 



(25) 

(26) 

(27) 
(28) 
(29) 

(30) 



(31) 
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with xe defined in Eq. (H). We find: 



1 dT 

To dxb 



2s 



Li2 



2Q 



- log(s + Q) log 



+ log 



1-s+Qj 
s-b + Q 



Li2 



2Q 



s-b + Q 



2s{l-P) 



+ 



1 



log(6) log 



+ (1-6) log 



fs-b + Q' 
l-s+Q\ 



+ (3g + .-5)lo, 



s + Q' 
Vb , 



w 



Q 



6 



{w — b){s — b) 



\og{b) 



2 logy 

2 



w 
1 



+ 



(5(1 - Xfe) 
' 1 + b' 



(1 - G, 



1 + 



2w 



+ 2sJxl-(3^ 



'2I 



I- Xb 



Go VI - 2sxb + b 



wGq 

{1-Xb)-1 
1 + b 



(32) 



2w 



{1-Xb)-1 



In Eqs. ( p4|) and (|3^) we note the presence of the so-called 'plus prescription' 1/(1 — X6) + . 
Such a term arises after one integrates over the phase space for real-gluon radiation. For 
a massive b quark, where Xb,min 7^ 0, one makes use of the following expansion: 



{i-XbY^^ 

with the plus prescription defined as: 



{Xb 3^6, mill) 1 e-z-i \ I 1 I /o/' ^ 

= -7^^' - + (TT^ + ^^^^ 



(33) 



1 



{g{xb)), h{xb)dxb 



'^b.min 



g{xb)[h{xb) - h{l)]dxb. 



(34) 



We observe that in Eq. (30) the quark mass regulates the collinear divergence, there- 
fore, unlike the result with a massless b quark (p4[), Eq. (30) does not contain any 
dependence on the dimensional- regularization quantities e or np. 

For the sake of checking the initial condition of the 6-quark perturbative fragmenta- 
tion function (p!oD , we need to rewrite Eq. (30) neglecting powers of mb/rrit. 

We find: 

1 dV , asifJ') 



Fo dxb 



S{l-Xb) + 



with 

Aiixb) 



1 -\- x'l 3 , , 



2tt 



-Ai{xb), 



(35) 



log — 
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l+x? Axb , Aw{l-w) Xbil-Xb) 
+ 2-- —\og[{l-w}Xb\--- 7— + 



[1-Xb)+ {1 - Xb)+ l + 2w l-{l-w)xb 

+ 5{l-Xb 



4Li2(l -w) + 2 log w log(l - w) - ^ - log(l - w) 

6 I + 2w 



2w 

log W — 4: 



1 — w 



(36) 



We note that in Eq. (|36|) the strong coupling constant multiplies the large logarithm 
log{mf/ml), which spoils the reliability of the fixed-order calculation and makes the 
approach of perturbative fragmentation mandatory. Such large logarithms are absent 
in the total NLO width, which can be checked to agree with the result in the literature. 



once one accounts for the mass of the b quark |36|. 

If fip is of the order of rrib, one can express the perturbative fragmentation function 
D{xb, fiF,^b) as a power expansion in as- 

D{xb, i^F, rub) = S^\xb) + ^^4^S^\xb, i^F, rub) + 0{al) (37) 

Inserting Eqs. (||), (|35D and (^) in Eq. (^), and solving for and one finds: 

S"\xb) = 6{1-Xb), (38) 
S'^\xb,^F,mb) = Ai{xb) - Ai{xb) ■ (39) 

Comparing then Ai{xb) and Ai{xb), it is straightforward getting Eq. (p^Ol). It should be 
noted that, although Ai{xb) and Ai{xb) separately depend on m^y via the ratio w, their 
difference and hence the initial condition for the perturbative fragmentation function do 
not, which is essential for it to be process independent. 

B Coefficient function in Mellin space 

For the sake of completeness, we wish to present the result for the Mellin transform of 
the MS coefficient function (Bl): 



MS 



fAr(mt,mvi/,/i,/iF) = / dz '^—^ {z,mt,mw, fi, fJ'F)- (40) 

1 "'0 az 

Given T{x) the Euler Gamma function, we define the poligamma functions: 



and the combinations 



Si{N) = MN+l)-M^), (42) 
S,{N) = -^i(iV+l) + ^i(l). (43) 
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The basic, non-trivial A^-space transforms of the terms which appear in Eq. (H) are given 
by: 



dz- — 

(1-^)+ 

dz-^z^-' 



dz 



1 - z 



1-z) 



logfl - z] 



1 z(l-z) 
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1 — (1 — w)z 



-S^{N-1), 

^[SI{N-1) + S2{N-1) 

2F^{l,N + l,N + 2,l-w) 
N+1 

2Fi{l,N + 2,N + 3,l-w) 
N + 2 



(44) 
(45) 
(46) 



(47) 



We shall then get: 
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l + 2w 
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+ Sl{N -l) + Sl{N + l) 



+ S2{N -l) + S2{N +l) + 2[l-2 log(l - w)]Si{N) 
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2 logwlog(l — w) — 2z — logfl — w) 



l + 2w 



2w 
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log w 



+ 4Li2(l-w7)-6- — |. 

One can show that, for N —>■ oo: 

Si{N) ~ ^o(A^) ~ log AT. 



(4J 



(49) 



Similar large- behaviour is also shown by the Mellin transform of the initial condition 
of the perturbative fragmentation function (p!OD, whose most-singular term at large Xb 
has a A^-space counterpart analogous to Eq. (|^). 
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